Orbital structure of the GJ876 extrasolar planetary system, based on the latest Keck and HARPS radial velocity data

We use full available array of radial velocity data, including recently published HARPS and Keck observatory sets, to characterize the orbital configuration of the planetary system orbiting GJ876. First, we propose and describe in detail a fast method to fit perturbed orbital configuration, based on the integration of the sensitivity equations inferred by the equations of the original $N$-body problem. Further, we find that it is unsatisfactory to treat the available radial velocity data for GJ876 in the traditional white noise model, because the actual noise appears autocorrelated (and demonstrates non-white frequency spectrum). The time scale of this correlation is about a few days, and the contribution of the correlated noise is about 2 m/s (i.e., similar to the level of internal errors in the Keck data). We propose a variation of the maximum-likelihood algorithm to estimate the orbital configuration of the system, taking into account the red noise effects. We show, in particular, that the non-zero orbital eccentricity of the innermost planet \emph{d}, obtained in previous studies, is likely a result of misinterpreted red noise in the data. In addition to offsets in some orbital parameters, the red noise also makes the fit uncertainties systematically underestimated (while they are treated in the traditional white noise model). Also, we show that the orbital eccentricity of the outermost planet is actually ill-determined, although bounded by $\sim 0.2$. Finally, we investigate possible orbital non-coplanarity of the system, and limit the mutual inclination between the planets \emph{b} and \emph{c} orbits by $5^\circ-15^\circ$, depending on the angular position of the mutual orbital nodes.Comment: 36 pages, 11 figures, 3 tables; Accepted to Celestial Mechanics and Dynamical Astronom

Resonances of low orders in the planetary system of HD37124

The full set of published radial velocity data (52 measurements from Keck + 58 ones from ELODIE + 17 ones from CORALIE) for the star HD37124 is analysed. Two families of dynamically stable high-eccentricity orbital solutions for the planetary system are found. In the first one, the outer planets c and d are trapped in the 2/1 mean-motion resonance. The second family of solutions corresponds to the 5/2 mean-motion resonance between these planets. In both families, the planets are locked in (or close to) an apsidal corotation resonance. In the case of the 2/1 MMR, it is an asymmetric apsidal corotation (with the difference between the longitudes of periastra $\Delta\omega\sim 60^\circ$), whereas in the case of the 5/2 MMR it is a symmetric antialigned one ($\Delta\omega = 180^\circ$). It remains also possible that the two outer planets are not trapped in an orbital resonance. Then their orbital eccentricities should be relatively small (less than, say, 0.15) and the ratio of their orbital periods is unlikely to exceed $2.3-2.5$.Comment: 28 pages, 10 figures, 3 tables; Accepted to Celestial Mechanics and Dynamical Astronom

On the dynamics of Extrasolar Planetary Systems under dissipation. Migration of planets

We study the dynamics of planetary systems with two planets moving in the same plane, when frictional forces act on the two planets, in addition to the gravitational forces. The model of the general three-body problem is used. Different laws of friction are considered. The topology of the phase space is essential in understanding the evolution of the system. The topology is determined by the families of stable and unstable periodic orbits, both symmetric and non symmetric. It is along the stable families, or close to them, that the planets migrate when dissipative forces act. At the critical points where the stability along the family changes, there is a bifurcation of a new family of stable periodic orbits and the migration process changes route and follows the new stable family up to large eccentricities or to a chaotic region. We consider both resonant and non resonant planetary systems. The 2/1, 3/1 and 3/2 resonances are studied. The migration to larger or smaller eccentricities depends on the particular law of friction. Also, in some cases the semimajor axes increase and in other cases they are stabilized. For particular laws of friction and for special values of the parameters of the frictional forces, it is possible to have partially stationary solutions, where the eccentricities and the semimajor axes are fixed.Comment: Accepted in Celestial Mechanics and Dynamical Astronom

Massive Search for Spot- A nd Facula-Crossing Events in 1598 Exoplanetary Transit Light Curves

We developed a dedicated statistical test for a massive detection of spot- A nd facula-crossing anomalies in multiple exoplanetary transit light curves, based on the frequentist p-value thresholding. This test was used to augment our algorithmic pipeline for transit light curves analysis. It was applied to 1598 amateur and professional transit observations of 26 targets being monitored in the EXPANSION project. We detected 109 statistically significant candidate events revealing a roughly 2 : 1 asymmetry in favor of spots-crossings over faculae-crossings. Although some candidate anomalies likely appear non-physical and originate from systematic errors, such asymmetry between negative and positive events should indicate a physical difference between the frequency of star spots and faculae. Detected spot-crossing events also reveal positive correlation between their amplitude and width, possibly due to spot size correlation. However, the frequency of all detectable crossing events appears just about a few per cent, so they cannot explain excessive transit timing noise observed for several targets